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Abstract 

We show that from the high order gravity theory it is possible to 
produce, in the linearized approch, particles which can be seen like massive 
modes of gravitational waves (GWs). The presence of the mass generates 
a longitudinal force in addition of the transverse one which is proper of 
the massless gravitational waves and the response an interferometer to the 
effect is computed. This could be, in principle, important to discriminate 
among the gravity theories. The presence of the mass could also have 
important applications in cosmology because the fact that gravitational 
waves can have mass could give a contribution to the dark matter of the 
Universe. 

PACS numbers: 04.80.Nn, 04.30.Nk, 04.50. +h 



1 Introduction 

The data analysis of interferometric GWs detectors has recently started (for the 
current status of GWs interferometers see [IllSISllHISlEllTllH]) and the scientific 
community hopes in a first direct detection of GWs in next years. The results of 
these detectors will have a fundamental impact on astrophysics and gravitation 
physics. There will be lots of experimental data to be analyzed, and theorists 
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will be forced to interact with lots of experiments and data analysts to extract 
the physics from the data stream. 

Detectors for GWs will be important for a better knowledge of the Universe 
and also to confirm or ruling out the physical consistency of General Relativity 
or of any other theory of gravitation [9l [lOl [TU [121 SSI [I^- This is because, in 
the context of Extended Theories of Gravity, some differences between General 
Relativity and the others theories can be pointed out starting by the linearized 
theory of gravity P [TOl [121 [14] . 

In this paper the production and the potential detection with interferometers 
of a hypotetical massive component of gravitational radiation which arises from 
the theory of gravity, which was the first and simplest high order gravity 
theory proposed [15], is shown. 

In the second Section of this paper it is shown that a massive mode of 
gravitational radiation arises from the high order action [15] 



Equation |(T]) is a particular choice with respect the well known canonical 
one of general relativity (the Einstein - Hilbert action [ElllT]) which is 



where R is the Ricci scalar curvature. We empahsize that the presence of 
the mass could also have important applications in cosmology because the fact 
that gravitational waves can have mass could give a contribution to the dark 
matter of the Universe. We also recall that an alternative way to resolve the 
dark matter and dark energy problems using high order gravity is shown in ref. 



In Section three it is shown that the massive component generates a longi- 
tudinal force in addition of the transverse one which is proper of the massless 
case. 

After this, in Section four, the potential interferometric detection of this 
longitudinal component is analyzed and the response of an interferometer is 
computed. This could be, in principle, important to discriminate among several 
gravity theories which are today considered. 

2 The production of a massive mode of gravita- 
tional radiation in the theory of gravity 

If the gravitational Lagrangian is nonlinear in the curvature invariants the Ein- 
stein field equations has an order higher than second [9l [121 ES] • For this reason 
such theories are often called higher-order gravitational theories. This is exactly 
the case of the action ([J). 




(1) 




(2) 
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By varying this action with respect to g^i, (see refs. \12\ [T3] for a parallel 
computation) the field equations are obtained (note that in this paper we work 
with G — 1, c — 1 and h = 1): 

_ -AttG f , rpjm) _ 1 ^d2 I 

(3) 

with associed a Klein - Gordon equation for the Ricci curvature scalar 

aR^m^{R + 8nGT), (4) 

where □ is the d' Alembertian operator and the mass m has been introduced 

1 
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for dimensional motivations: = — t^, thus a has to be negative [15] 



In the above equations T^™"* is the ordinary stress-energy tensor of the matter 
and G is a dimensional, strictly positive, constant [HI [El SI]. The Newton 
constant is replaced by the effective coupling 

which is different from G. General relativity is obtained when a = 0. 

To study gravitational waves the linearized theory in vacuum (T^T^ = 0) has 
to be analyzed, with a little perturbation of the background, which is assumed 
given by the Minkowskian background. In this case the Ricci scalar is assumed 
slowly varying near zero: i? ~ + SR = hn. 

Putting 

to first order in hf^^, , calling R^ppa , R^u and R the linearized quantity which 
correspond to Rfiupa , R^lv and R, the linearized field equations are obtained 
fHfllfllflTl: 



(7) 



Riiy - ^ritiu = d^dyR + i-jp^niiR 
R^ivpa and eqs. lO are invariants for gauge transformations [T2t [T3] 

hn h'j^ ~ hji; 

then 



(8) 



hfj^i, = hfj^^ - ^rj^i, + ri^i,hR (9) 
can be defined, and, considering the transform for the parameter 
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□e. = d^h^^, (10) 
a gauge parallel to the Lorenz one of electromagnetic waves can be choosen: 

d^h^. = 0. (11) 
In this way field equations read like 

Dh^, = (12) 

OhB. = m'^hR (13) 
Solutions of eqs. ifT2|) and (fT3)) are plan waves: 



hf,i, = ^^L/("^) expiip^Xa) + c.c. (14) 
h-R = a(~P*) exp{iq°'xa) + c.c. (15) 



where 



{uj,^) u=p=\'p'\ 



(16) 



In eqs. (|T2|) and lfT4| the equation and the solution for the waves like in 
standard general relativity [I6l[l7] have been obtained, but eqs. (fT3| and lfT5|) 
are respectively the equation and the solution for the massive mode (see also 
\12\ 113)) arising from the Starobinsky's high order gravity theory. 

The fact that the dispersion law for the modes of the massive field Hr is not 
Hnear has to be emphatized. The velocity of every tensorial mode h^j,i, is the 
Hght speed c, but the dispersion law (the second of eq. (fT6|) ) for the modes of 
Kr is that of a massive field which can be discussed like a wave-packet O US] . 
Also, the group-velocity of a wave-packet of Hr centered in "p* is 

= (17) 

which is exactly the velocity of a massive particle with mass m and momen- 
tum 7^. 

From the second of eqs. (fT6| and eq. I|17p it is simple to obtain: 



VG = • (18) 

Then, wanting a constant speed of our wave-packet, it has to be [13] 

= \/(1-^^gV- (19) 
The relation (fT9|) is shown in fig. 1 for a value vg = 0.9. 
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Figure 1: the mass-frequency relation for a massive graviiatonal wave arising 
from the high order gravity theory and propagating with a speed of 0.9c : 
for the mass it is IHz = 10~^^eV 



Now the analisys can remain in the Lorenz gauge with trasformations of the 
type De^ = 0; this gauge gives a condition of transversahty for the tensorial 
part of the field: k^A^i^ = 0, but does not give the transversahty for the total 
field hfj^^. From eq. ^ it is 

hf^i, = hf^^ - -^Jy^L/ + Tj^yhR. (20) 
At this point, if being in the massless case [17], it could been put 

□e'' = 

(21) 

a^e^ = -'i + hn, 

which gives the total transversahty of the field. But in the massive case this 
is impossible. In fact, applying the Dalembertian operator to the second of eqs. 
l(2T|) and using the field equations lfT2|) and (fTSl it results 

□e^ = m'^hR, (22) 

which is in contrast with the first of eqs. ((2T1) . In the same way it is possible 
to show that it does not exist any linear relation between the field h^y and Hr. 
Thus a gauge in wich h^i, is purely spatial cannot be chosen (i.e. it cannot be 
put /ipo — 0, see eq. I|20p ) . But the traceless condition to the field h^y can be 
put: 

□e'' = 

(23) 

These equations imply 

^^'h^,, = 0. (24) 
To save the conditions d^.h^'^ and h = transformations like 

□e^ = 

(25) 

= 

can be used and, taking in the z direction, a gauge in which only An, 
A22, and A12 ~ A21 are different to zero can be chosen. The condition h = 
gives All ~ —A22. Now, putting these equations in eq. I|20p it results 

V(t, z) = A+{t - z)el+J + yl^(t - z)e|,^) + hait - VGz)ri^^. (26) 

The term [t — z)elt}) + A^ (t— z)e"f^J describes the two standard polariza- 
tions of gravitational waves which arise from General Relativity, while the term 
hR{t — VGz)r]fj_i, is the massive polarization arising from the theory. 
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3 The presence of a longitudinal force 



The analysis of the two standard polarization is well known in the literature 
[ElIlT]- For a the pure polarization arising by the theory eq. (|26)) can be 
rewritten as 

hf,^{t - vgz) = hfiit - VGz)r]f^^ (27) 
and the corrispondent line element is the conformally flat one 

ds^ = [1 + hnit - VGz)]{-dt^ + dz^ + dx^ + dy^). (28) 

But, in a laboratory environment on Earth, the coordinate system in which 
the space-time is locally flat is typically used and the distance between any 
two points is given simply by the difference in their coordinates in the sense of 
Newtonian physics [H HSl [lH [17] . This frame is the proper reference frame of 
a local observer, located for example in the position of the beam splitter of an 
interferometer. In this frame gravitational waves manifest themself by exerting 
tidal forces on the masses (the mirror and the beam-splitter in the case of an 
interferometer) . A detailed analysis of the frame of the local observer is given 
in ref. [17], sect. 13.6. Here only the more important features of this coordinate 
system are recalled: 

the time coordinate xq is the proper time of the observer 0; 

spatial axes are centered in 0; 

in the special case of zero acceleration and zero rotation the spatial coor- 
dinates the proper distances along the axes and the frame of the local 
observer reduces to a local Lorentz frame: in this case the Hne element reads 

m 

ds^ = -{dx^f + S.jdx'dx^ + 0{\x^^)dx°'dx'^. (29) 

The effect of the gravitational wave on test masses is described by the equa- 
tion 

? = -Rlkox\ (30) 

which is the equation for geodesic deviation in this frame. 
_ Thus, to study the effect of the massive gravitational wave on test masses, 
Rqi^q has to be computed in the proper reference frame of the local observer. 
But, because the linearized Riemann tensor R^i,pa is invariant under gauge 
transformations [l2l SHI [17] , it can be directly computed from eq. l(27j) . 
From [TT] it is: 

Rfiupa ~ -^{dfj,d(3hau + dudah^fj — dadph^i, — dpdvhap}, (31) 
that, in the case eq. (|27|) . begins 

^070 = \{d''dahnm^ + d^d^hnd^ - d'^d^liRrioa ~ dodohnS^}; (32) 
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the different elements are (only the non zero ones will be written): 



dod^hnSQ 



dfhji for a = 7 = 



—dzdthfi for q; = 3;7 = 
dfhfi for a = 7 = 

dtdzhu for q; = 0;7 = 3 



(33) 
(34) 



Now, putting these results in eq. ([32|) it results: 



-dfhR 


for 


a = 7 = 


^ 


dlhR 


for 


a — 7 = 


3 


-dtdzhR 


for 


a = 0; 7 = 


= 3 


dzdthR 


for 


a = 3;7 = 


= 


-dlhR 


for a 


= 7 ■ 





-^010 - ~^^f 



-^01 



-^hR 



-^030 - \^^R- 



But, putting the field equation (fT3l) in the third of eqs. l(37|) it is 



r3 

^030 



1 2, 

2"^- riR, 



which shows that the field is not transversal. 
Infact, using eq. ([3Q| it results 



1. 



-J^rx, 



and 



z = ^Ri^ ~ vgz)z. 



(35) 



(36) 



(37) 



(38) 

(39) 
(40) 

(41) 



Then the effect of the mass is the generation of a longitudinal force (in 
addition to the transverse one). Note that in the Hmit to ^ the longitudinal 
force vanishes. 
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4 The interferometer's response to the longitudi- 
nal component 

Before starting the analysis it has to be discussed if there are fenomenogical 
limitations to the mass of the wave OIIS]. Treating like a classical wave, 
that acts coherently with the interferometer, it has to be m ^ 1/L , where 
L = 3 kilometers in the case of Virgo and L = 4 kilometers in the case of LIGO. 
Thus it has to be approximately m < 10~^eV. However there is a stronger 
limitation coming from the fact that the massive wave needs a frequency which 
falls in the frequency-range for earth based gravitational antennas that is the 
interval lOHz < f < \QKHz [HEllSlllllllEllIllH]. For a massive gravitational 
wave, from the second of eqs. (fT6| it is: 



27r/ = w = yj-m? (42) 
were p is the momentum [13] . Thus it needs 

QeV <m< IQ-^^eV. (43) 

For these light particles their effect can be still discussed as a coherent grav- 
itational wave. For the discussion of this longitudinal effect we start directly 
from the gauge ([28|l . 



Eq. ([28| can be rewritten as 

where r is the proper time of the test masses. 

From eqs. l|28l) and (|44l) the geodesic equations of motion for test masses 
(i.e. the beam-splitter and the mirrors of the interferometer), can be obtained 



dr2 







(45) 

d^t _ 1 dt{l+hR) 

dP^ - 2 {l + hR^ 

£z_ _ 1 8.(l+feH) 



The first and the second of eqs. (|45|) can be immediately integrated obtaining 

Ci = const. (46) 



dx 
d7 



^ C2 = const. (47) 
or 
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In this way eq. (|44l) becomes 



If we assume that test masses are at rest initially we get Ci = C2 = 0. Thus 
we see that, even if the GW arrives at test masses, we do not have motion of 
test masses within the x — y plane in this gauge. We could understand this 
directly from eq. ((28l) because the absence of the x and of the y dependences in 
the metric impHes that test masses momentum in these directions (i.e. Ci and 
C2 respectively) is conserved. This results, for example, from the fact that in 
this case the x and y coordinates do not esplicitly enter in the Hamilton-Jacobi 
equation for a test mass in a gravitational field |16j . 

Now we will see that, in presence of the GW, we have motion of test masses 
in the z direction which is the direction of the propagating wave. An analysis 
of eqs. (|45|) shows that, to simplify equations, we can introduce the retarded 
and advanced time coordinates (w, w): 

u = t — VqZ 

(49) 

V =t + VQZ. 



From the third and the fourth of eqs. (|45l) we have 
d du dy [1 + hft (u)] 



dT dr (1 + hniu))^ 
This equation can be integrated obtaining 



0. (50) 



where a is an integration constant. From eqs. (|48| and (jSTj). we also get 

dv 13 



dr \ + hi- 

where B = and 



(52) 



T = /3u + 7, (53) 

where the integration constant 7 correspondes simply to the retarded time 
coordinate translation u. Thus, without loss of generality, we can put it equal 
to zero. Now let us see what is the meaning of the other integration constant 
/3. We can write the equation for z from eqs. fSTj) and l(52l) : 



dz 1 



10 



When it is ft-fl = (i.e. before the GW arrives at the test masses) eq. I|54p 
becomes 

But this is exactly the initial velocity of the test mass, so we have to choose 
P = 1 because we suppose that test masses are at rest initially. This also imply 
a = 1. 

To find the motion of a test mass in the z direction we see that from eq. 
l(53|) we have dr = du, while from eq. l(52]) we have dv = j^^- Because it is 
vg^ — we obtain 



i , dr , ^ , , 



which can be integrated as 



(57) 



where zq is the initial position of the test mass. Now the displacement of 
the test mass in the z direction can be written as 



(58) 



(59) 



Az^Z-Zo^-^ f'-^^^o-voAz hnH^ 

^ _ 1 rt^-a^o ^^du. 
We can also rewrite our results in function of the time coordinate t: 

x{t) = Xq 

T{t) = t - VGz{t), 

Calling I and L + / the unperturbed positions of the beam-splitter and of the 
mirror and using the third of eqs. I|59p the varying position of the beam-splitter 
and of the mirror are given by 

(60) 

-M(t) = L + /-^/!--^^+"^d(u) 

But we are interested in variations in the proper distance (time) of test 
masses, thus, in correspondence of eqs. ijGO]). using the fourth of eqs. I|59p we 
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get 

Then the total variation of the proper time is given by 



(61) 



1 rt-VG(L+l) hr,(u) 

A T{t) = TM{t) - TBs{t) = vgL-- , A ^)- (62) 



2 Jt^vGi 1 + ^«(") 

In this way, recaUing that in the used units the unperturbed proper distance 
(time)is T = L, the difference between the total variation of the proper time in 
presence and the total variation of the proper time in absence of the GW is 

6r{t) EE Ar{t) -L = -L{vg + I) - .J", M - (63) 

This quantity can be computed in the frequency domain, defining the Fourier 
transform of hn as 



hniuj) ^ / dt hR{t)cyij){iujt). (64) 
and using the translation and derivation Fourier theorems, obtaining 

6?{uj) = L(l - vl) exp[iujL{l + va)] + 2a>L(4-i)^ 

[ex.p[2iujL]{vG + lf{-2i + LuL{vG - 1) + 2L cxp[iujL{l + vg)] (65) 

{6ivG + 2iv% -ujL + uoLv%) + L{vg + l f{-2i + ljL{vg + l))]^ii- 
A "signal" can be also defined: 

Siu;) ^^ = il-vl) exp[*c.L(l + vg)] + 

[cxp[2iLuL] [vG + lf{-2i + ujL{vg - 1) + 2 cxp[iwL(l + vg)] (66) 

(6wg + 2wc; -uoL + ujLv%) + [vg + l)^(-2i + ujL{vg + ^))\hR- 
Then the function 

Tiioj) = (1 -i;^)cxp[iwL(l + WG)] + 2u>L(vl-iy 

[cxp,[2iujL\{vG + lf{-2i + ujL(vG - 1) + 2 cxp[iwL(l + wg)] (67) 
{QivG + 2ivl, -ujL + ujLv%) + [vg + l)^(-2i + ujL{vg + 1))], 
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Figure 2: the absolute value of the longitudinal response function l(65|) of the 
Virgo interferometer (L = iKm) to a GW arising from the high order gravity 
theory and propagating with a speed of 0.1c (non relativistic case). 

is the response function of an arm of our interferometer located in the z-axis, 
due to the longitudinal component of the massive gravitational wave arising from 
the high order gravity theory and propagating in the same direction of the 
axis. 

For vg ^ 1 it is Ti(w) 0. 

In figures 2, 3 and 4 are shown the response functions l(67|) for an arm of the 
Virgo interferometer (L = 3Km) for vg = 0.1 (non-relativistic case), vg = 0.9 
(relativistic case) and vg = 0.999 (ultra-relativistic case). We see that in the 
non-relativistic case the signal is stronger as it could be expected (for m ^ 
we expectTi(a;) ^0). In figures 5, 6, and 7 the same response functions are 
shown for the Ligo interferometer (L — AKm). 

5 Conclusions 

We have shown that from the R^ high order gravity theory it is possible to pro- 
duce, in the linearized approch, particles which can be seen like massive modes 
of gravitational waves. The presence of the mass generates a longitudinal force 
in addition of the transverse one which is proper of the massless gravitational 
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Figure 3: the absolute value of the longitudinal response function l(65|) 
Virgo interferometer (L = 3Km) to a GW arising from the high order 
theory and propagating with a speed of 0.9 (relativistic case). 
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Figure 4: the absolute value of the longitudinal response function l(65|) of the 
Virgo interferometer (L = SKm) to a GW arising from the high order gravity 
theory and propagating with a speed of 0.999 (ultra relativistic case). 
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Figure 5: the absolute value of the longitudinal response function ([65]) of the 
LIGO interferometer {L ~ AKm) to a GW arising from the high order 
gravity theory and propagating with a speed of 0.1c (non relativistic case). 
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Figure 6: the absolute value of the longitudinal response function ijGS]) of the 
LIGO interferometer {L = AKm) to a GW arising from the high order 
gravity theory and propagating with a speed of 0.9c (relativistic case). 
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Figure 7: the absolute value of the longitudinal response function l(65|) of the 
LIGO interferometer {L = AKm) to a GW arising from the high order 
gravity theory and propagating with a speed of 0.999c (ultra relativistic case). 
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waves and the response an interferometer to the effect has been computed. The 
presence of the mass could also have important applications in cosmology be- 
cause the fact that gravitational waves can have mass could give a contribution 
to the dark matter of the Universe. As a final remark, we recall that the po- 
tential detection of a longitudinal component of GWs could be, in principle, 
an useful tool to discriminate among several gravity theories which are today 
considered. 
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